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We study the local structure of zero mode wave functions of chiral matter fields in 



from local Higgsing in the 8-dimensional parent action of F-theory 7-branes. The 
solutions are found as expansions both in powers and derivatives of the magnetic 
fluxes. Yukawa couplings are given by an overlap integral of the three wave func- 
tions involved in the interaction and can be calculated analytically. We provide 
explicit expressions for these Yukawas to second order both in the flux and deriva- 
tive expansions and discuss the effect of higher order terms. We explicitly describe 
the dependence of the couplings on the U(l) charges of the relevant fields, appro- 
priately taking into account their normalization. A hierarchical Yukawa structure 
is naturally obtained. The application of our results to the understanding of the 
observed hierarchies of quarks and leptons is discussed. 



1 Introduction 



The hierarchical structure of fermion masses and mixings is one of the most remark- 
able properties of the Standard Model (SM). An outstanding challenge of string theory 
compactifications is to obtain models with the massless spectrum of the SM and repro- 
ducing naturally such hierarchical structure. In type IIB orientifold, as well as heterotic, 
compactifications the Yukawa couplings which govern fermion masses and mixings are in 
principle calculable, in the large compact volume limit, in terms of overlap integrals [1,2] 



Here ipi and <pH are internal wave functions associated to the fermions and Higgs fields 
respectively, taking values in the compact complex threefold X 3 . These wave functions 
are zero modes of higher dimensional fields in the compact internal space. The technical 
problem here is that in general we do not know how to compute the relevant wave functions 
for arbitrary curved spaces X 3 . Such a computation has been completely worked out for 
the relatively simple case of type IIB toroidal orientifolds [3] with constant U(l) fluxes 
(see also [4-6]). In this case with a flat geometry the equations of motion can be fully 
solved to obtain the wave functions which turn out to have a neat expression in terms 
of Jacobi ^-functions. It was found that in the simplest models only one generation of 
quarks and leptons acquires a non-trivial Yukawa coupling, which is a good starting point 
to reproduce the observed hierarchies [3] . Having explicit solutions for the wave functions 
is also useful to study other physical properties of the compactifications such as the effect 
of closed string fluxes and warping [7]. 

Clearly, it would be interesting to obtain wave functions and Yukawa couplings in more 
complicated curved geometries and for non-constant fluxes. An obvious obstruction is that 
determining the wave functions seems to require a knowledge of the global geometry of 
the compact X 3 manifold. In fact, the problem may be more tractable within the context 
of a bottom-up approach as advocated in [8] (see also [9-12]). The idea is that in order 
to extract the relevant physics of a SM compactification it is enough to have a local 
description of the geometry of the branes in which the SM fields reside. This is the case 
for example of models derived from D3-branes at singularities [8,9,11] in which the SM 
physics only depends on the local geometry around the singularity. This type of structure 
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Figure 1: Intersecting matter curves 



is also characteristic of local configurations of D7-branes wrapping intersecting 4-cycles 
inside X 3 . F-theory [13] is the natural non-perturbative extension of these local 7-brane 
configurations. In the last year local F-theory GUT constructions have been proposed 
[14-16] as a particularly attractive class of bottom-up configurations with a number of 
phenomenological virtues (see [17-31], as well as [32-43] for other recent developments). 
In this scheme the Yukawa couplings arise again as overlap integrals now of the form 



in which S is the compact complex twofold wrapped by the GUT F-theory 7-brane. The 
quark and lepton multiplets of the SM reside at matter Riemann curves £/ inside S, 
which correspond geometrically to the intersection of S with the world-volume of other 
U(l) 7-branes. Yukawa couplings come from the triple overlap of these matter curves 
involving quarks, leptons and Higgs fields (see figure 1). In order to compute the Yukawa 
coupling (1.2) we need again the internal wave functions. However, in this case given 
the local geometry of the coupling it would be enough to have a knowledge of the wave 
functions close to the intersection point. It was pointed out in [15] that one can determine 
the profile of these wave functions close to the intersection point in terms of a certain 
quasi-topological theory in D=8. The equations of motion of that theory have solutions 
corresponding to hypermultiplet zero modes localized along the matter curves with a 
Gaussian profile. One finds [19] that to leading order, for a compactification having three 
generations, only one of them gets a non-trivial Yukawa, in analogy with the results in [3] . 




(1.2) 



However the distortion of the wave functions, due to the presence of U(l) gauge fluxes, 
could be the natural source of the observed hierarchy of masses and mixings of quarks 
and leptons [19]. 

In this article we make a systematic study of the solutions of the differential equations 
of motion of the quasi-topological D=8 theory of [15]. The zero mode solutions give the 
local wave functions corresponding to the massless particles residing at the matter curves 
of F-theory unification models. We make an expansion both in powers and derivatives 
of the U{1) fluxes and explicitly solve the differential equations. Equipped with these 
wave functions we compute the Yukawa couplings from the overlap integral of the three 
wave functions involved in the couplings. These integrals may be calculated analytically 
and we provide explicit expressions for these Yukawa couplings up to fourth order in the 
flux and derivative expansions. As suggested in [19], a hierarchy of masses for fermions 
naturally appears. We also study the application of our results to the understanding of the 
observed hierarchy of masses and mixings in the SM. We find good qualitative agreement 
with experiment for reasonable ranges of the flux parameters. 

The organization of the rest of this article is as follows. In the next section we provide 
a brief review of the aspects of F-theory models that concern our discussion. In chapter 3 
we study the wave functions of the zero modes which are solutions of the quasi-topological 
D=8 field theory equations. We consider both constant and varying fluxes in a general 
setting of three intersecting matter curves. The details of the solutions are given in 
appendix A. In chapter 4 we address the explicit computation of the Yukawa couplings 
by evaluating the overlap integral of the three relevant wave functions. Based on the 
leading terms in the Yukawa couplings provided in appendix B, we describe the general 
structure of the flux-induced corrections and their contribution to the Yukawa matrices. 
In chapter 5 we apply the previous results to the analysis of the fermion mass spectra of 
SU(5) GUT's with non-vanishing hypercharge flux breaking the theory down to the SM. 
We show that reasonable agreement with observed mass hierarchies and mixing may be 
obtained for appropriate flux parameters. Chapter 6 is devoted to some final comments 
and discussion. 

The computation of Yukawa couplings in analogous settings has been more recently 
analyzed in [49,50], see the note added at the end of this paper. 



2 Review of F-theory unification 



The purpose of this section is to give a short overview of the F-theory formalism developed 
in [15] (see also [14,17]). 

In the F-theory setup, the D=4 supersymmetric gauge theory descends from 7-branes 
wrapping a compact surface S of complex codimension one in the threefold base of an 
elliptically-fibered Calabi-Yau fourfold. The gauge group G s on the 7-branes depends on 
the singularity type of the elliptic fiber. In turn Gs can be broken by a vev in a subgroup 
Hs C Gs- We consider Hs = U(l), typical examples being the hypercharge in SU(5) or 
U(1)b-l m an 5*0(10) GUT. The U(l) background breaks the gauge group and gives rise 
to matter charged under the commutant of H s in Gs- We assume that S is a del Pezzo 
surface so that gravity decouples from the gauge theory [15]. 

The singularity type of the elliptic fiber can be enhanced to group G^ along a curve 
£ C S of complex codimension two on the threefold base. This curve appears at the 
intersection of S and another surface S'. On the 7-branes wrapping S' there is a gauge 
theory with group Gs> which decouples when S' is non-compact. Based on the knowledge 
of intersecting D-branes, one expects additional degrees of freedom due to open strings 
stretching between the 7-branes wound on S and S'. The extra fields localized on the 
matter curve £ must be charged under Gs x Gs>- Indeed this is the picture that arises in 
F-theory [15,44]. 

We will now review the basic facts about the charged fields originating in the surface 
S and in the matter curve £. Our discussion is brief and follows mostly [15]. 

2.1 Bulk fields 

The effective physics of the 7-branes wrapping S is described by D=8 twisted super 
Yang-Mills on R 3 ' 1 x S [15]. The supersymmetric multiplets include the gauge field, 
plus a complex scalar tp and fermions (i],ifj,x) in the adjoint. After twisting the scalar 
and fermions become forms on S. Using local coordinates (zi,z 2 ) for S the results are 
summarized by 

A = A^dx 1 * + A m dz m + Afndz m ; y? = (p 12 dz 1 A dz 2 

if>a = ikai dz 1 + i\) a2 - dz 2 ; Xa = Xan dz 1 A dz 2 (2.1) 



Notice that ip is a (0,1) form whereas ip and \ are (2,0) forms. The remaining fermion r) a 
is a (0,0) form. The subscript a, which corresponds to left handed fermions in R 3,1 , will 
be dropped hereafter. The D=4, Af=l theory has gauge multiplet (A^,rj), together with 
chiral multiplets (A m ,ip m ) and (y?i2, X12)) phis their complex conjugates. 

The D=8 effective action found in [15] can be integrated over the compact surface S 
to obtain the dynamics of the D=4 multiplets. In computing couplings of the charged 
fields the most interesting term will be the superpotential 

W = Mt I Tr {Ff 2) A $) = [ Tr (3 A A $) + [ Tr (A A A A $) (2.2) 
Js Js Js 

where M* is the mass scale characteristic of the supergravity limit of F-theory. Here A 

and $ are chiral superfields with components 

A m = 4 + V / 2^m + '-- 

$ 12 = tp 12 + V29 X i2 + ■ ■ ■ (2.3) 

where • • • involves auxiliary fields. Only the (0,2) component of the superstrength appears 
in (2.3). 

The equations of motion derived from the D=8 effective action are the starting point 
to discuss the zero modes. The part of the action bilinear in fermions, without kinetic 
terms, is found to be [15] 

I F = I d 4 x Tr ( X A d A ij + 2iV2u A d A r] A ip + -if> A [<p, ip] + V2r)[<p, x] + h.c.) (2.4) 

where ui is the fundamental form of S. Taking variations with respect to 77, ip and x 
respectively gives the equations of motion 

uAd A ^ + ^,x] = (2-5) 

d A X - 2iV2u Adrj- [ip, ip] = (2.6) 

8 A il>-V2[<p,rj\ = (2.7) 

For the bosonic fields it is found that the field strength F s must have vanishing (2,0) and 
(0,2) components and verify the BPS condition 

ujAF s + l -[ip^\ = (2.8) 



Finally, the complex scalar must satisfy the holomorphicity condition 8aP = 0. 

To determine the charged massless multiplets in D=4 it is necessary to specify the 
background for the adjoint scalar ip and the gauge field. When (ip) = 0, the equations 
of motion imply that the number of zero modes of ip and x are counted by topological 
invariants that depend both on S and the gauge bundle of the background [15]. 

2.2 Fields at intersections 

We now want to discuss the degrees of freedom localized on a matter curve E occurring 
at the intersection of surfaces S and S'. As explained in [15], to preserve Af=l supersym- 
metry in D=A, the theory on R 3 ' 1 x E must be D=6 twisted super Yang-Mills. The D=6 
twisted supermultiplet, which includes two complex scalars and a Weyl spinor, decom- 
poses into D=A chiral multiplets (a, A) and (a c , A c ), plus CPT conjugates. The number 
of zero modes is given by topological invariants depending on E and the gauge bundle on 
a background in G^. 

There is a very nice intuitive way of understanding the matter localized at E. The 
idea, originally given in [44] and expanded in [15], is to start from the D=8 theory on S 
with gauge group Gs and then turn on a background for the adjoint scalar given by 

(<p) = m 2 z 1 Q 1 (2.9) 

where z\ is a complex coordinate on S, and Q\ is a Z7 (l)i generator in the Cartan sub- 
algebra of Gj> To streamline notation, p = ip 12 . We have explicitly introduced a mass 
parameter m so that p has the standard dimensions. The basic idea is that in presence 
of (p) the D=8 fields have zero modes localized at z\ — that are naturally associated 
to the fields at the intersection. 

When Z\ — the gauge group is unbroken, but when z\ ^ the group is broken to 
GsxC/(l)i, with G s being the group whose generators commute with Q x . The locus z\ — 
defines the curve E. Thus, on E the singularity enhances from Gs to Gs D Gs x U(l)i. 
The breaking of the gauge group is explained by the deformation of the singularity type 
from Gs to Gs due to the background in the Cartan subalgebra [44]. 

For ordinary D7-branes the adjoint scalar corresponds to degrees of freedom in the 
transverse direction and a non-zero vev means that some branes are separated and the 



gauge group is broken. For instance, if there are (K + 1) D7-branes to begin and one is 
moved away, the original SU(K + 1) is broken to SU(K) x U(l). Furthermore, the open 
strings stretching between the two stacks of D7-branes give rise to massless bifundamen- 
tals (K, —1) + (K, 1) localized at the intersection. For F-theory seven-branes wrapping 
a surface 5 one then expects (</?) to break the original gauge group to some Gs- More- 
over, there will be massless 'bifundamentals' descending from the adjoint of G s which 
decomposes as a direct sum of irreducible representations (R,q±) under Gs x U(l)i. 

Several examples of singularity resolution were worked out in [44] and more recently 
in [15,16,34]. For illustration let us consider the case G s = E 6 and G s = SO (10) that 
will be of interest later on. Under SO (10) x U(l) the E e adjoint decomposes as 

78 = (45, 0) + (1, 0) + (16, 1) + (16, -1) (2.10) 

Therefore, there will be chiral multiplets transforming as 16 and 16 of SO (10). To see 
how 5*0(10) is enhanced to E e on E it is convenient to represent the Cartan generators 
as vectors \Qi) so that \(p) oc z±\Qi) corresponds to the adjoint vev [44]. The simple roots 
are elements (vj\ of the dual space. Those roots with (vj\ip) = remain as 50(10) roots 
while those with (vj\tp) oc z\ become the weights of the 16 and 16. 

The resolution of the singularity by the adjoint vev can be figured out as explained 
in [44]. The generic E e singularity can be cast as [15,45] 

y 2 = x 3 + ^z 4 + e 2 xz 2 + e 5 xz + e 6 z 2 + e 8 x + e 9 z + e 12 ( 2 -H) 

where the are functions that depend on the adjoint vev. More precisely, in [45] they 
are given in terms of an arbitrary vector (t±, • • • ,te) in the Eq Cartan subalgebra. In 
our notation, in the E 6 D 5O(10) x U(l)i example, t\ = z\ while other t £ 's vanish. By 
choosing t\ = —St, and computing the according to the formulas of [45], we obtain the 
deformation 

y 2 = x 3 + V - 3t 2 xz 2 - 12t 5 xz - 6t e z 2 - 12t s x - 16t 9 z - 12t 12 (2.12) 

This is the same result found in [44], for a different though equivalent choice of vev vector. 
It can be shown that for t ^ there is an 5O(10) singularity. 

So far we have just reviewed how the gauge group on the curve E is enhanced. We 
now want to discuss how the matter localized on E arises from zero modes of the D=8 



bulk fields. It is enough to look at fermions because the scalars follow by supersymmetry. 
We then want to solve the D=8 equations of motion for the twisted fermions when f has 
the vev linear in zi, and there is no gauge background. The solutions that are localized 
at z\ — can be interpreted as the fermions A and A c that come from the twisted super 
Yang-Mills on R 3 ' 1 x S. 

We start from the D=8 fermionic equations of motion (2.5-2.7). To show that there 
are localized solutions it suffices to work locally and assume that the fundamental form 
of S has the canonical Euclidean form 

uj = Udz 1 Adz x + dz 2 Adz 2 ) (2.13) 

Notice that the coordinate along E is z 2 whereas the transverse coordinate is z\. To look 
for localized solutions one can neglect derivatives in z 2 . The equations of motion reduce 
then to 

V2 dif) — mfziqifa — ; d\^2 — V2m 2 ziqii] = (2-14) 
<9i^i - m 2 z x q x x = ; hx~ m 2 z 1 q 1 ^- 1 = (2.15) 

where \ — Xi2- Here q 1 is the U(l)i charge of the fermions that belong to a representation 
(R, q\) oiGs><U{l)i. From the above equations we see that there are no localized solutions 
for 77 and ip2, an d indeed it is consistent to set i] = and ipo. — 0- On the other hand, the 
coupled system for x an d tpi has solution 

X = f{z 2 )e-* m2 W ; = ~f(z2) e"«» 2 l«l a (2.16) 

where f(z 2 ) is an arbitrary holomorphic function of the coordinate along S. We have 
set qi — e to take into account normalization of the charges. Clearly the zero modes are 
peaked around z\ = 0, with width in \z±\ 2 equal to 1/em 2 . The constant em 2 is expected 
to be of the order of the F-theory mass scale M 2 . 

The solutions localized at z\ — naturally correspond to the fermions A and A c that 
appear in the D=6 twisted theory. As argued in [15], the transformations along S of 
(■01, x) and (A, A c ) do agree. 



3 Zero modes at intersecting matter curves 

As we have reviewed, there are charged fields localized on a matter curve E where the 
singularity type is enhanced. In this section we want to study the situation in which there 
are three matter curves £/, I — a,b, c, occurring at the intersection of surfaces S and S'j. 
In turn the three matter curves intersect at a point. On each curve there is a gauge group 
Gs 7 that enhances to G p at the common point of intersection [15]. 

To obtain the wave functions of the fermionic zero modes living at intersections we 
follow again the approach of [15]. The strategy is to consider the fermionic equations of 
motion of the D=8 theory on S with a non-trivial background for the adjoint scalar tp 
that determines the curves. One then looks for solutions that are localized on a particular 
matter curve. 

In the previous section we have seen that the equations that can give rise to localized 
fermionic zero modes are given by 



We have set r] = because only fermions that appear in D=A chiral multiplets are 
expected to have localized modes on E. Notice then that equation (2.7) implies the 
additional constraint Ba^P = 0. The new ingredient now is a more general background for 
the adjoint scalar tp. Concretely, 



where the Zi are local coordinates, and tp = tp l2 . Each Qi is the generator of a U(l)i 
inside G p . The mi are some mass parameters expected to be related to the F-theory 
supergravity scale M*. In what follows we will take m\ = m 2 = m. 

As discussed in section 2.2, when (tp) oc z±, the adjoint vev resolves the singularity 
on the curve E a characterized by z\ = [44]. Now the more general adjoint background 
resolves the G p singularity where three curves intersect. When (zi,z 2 ) 7^ (0,0), the 
group is broken to Gs but at the intersection (-21,-22) = (0,0) the group is enhanced to 
G p . Furthermore, at the curves £/ the group is enhanced to G S/ D G s x U(l). For 
example, when (G p ,Gs) = (£7, 50(10)), the group is enhanced to E 6 at E a and £5 
defined respectively by the loci z\ — and z 2 = 0, whereas it is enhanced to 50(12) at 



uAd A rl> + -[{$), X ] = 





(tp) = ml zi Qi + ml z 2 Q2 



(3.2) 



S c defined by z\ + z 2 = 0. In the case (G p , Gs) = {E 8 , E 6 ), the group is enhanced to E 7 
at each curve. 

At each curve X/ there are matter fermions that correspond to open strings stretching 
between 7-branes wrapping S and S'j. The U(l)i charges of these fermions, denoted 
(qi,q 2 ), depend on the curve as shown in table 1. In this table we also indicate how the 
fermions transform under the gauge group in the examples Gs = E 6 and Gs = 5*0(10) in 
which the group G p is respectively E s and E 7 . For Gs = SU(5) the rank two higher G p 
can be either E 6 or SO{12). We have introduced parameters (ei, e 2 ) to take into account 
normalization of the charges. 



curve 


(91, 92) 


locus 


_E 6 


SO(10) 


S77(5) 


SU(5) 


y 


(ei,0) 


zi = 


27 


16 


10 


10 




(0,e 2 ) 


z 2 = 


27 


16 


10 


5 




(-ei,-e 2 ) 


-21 + -22 = 


27 


10 


5 


5 



Table 1: Curves and charges 



3.1 Zero modes in the absence of fluxes 

We will first solve the zero mode equations without turning on a background gauge field 
but with scalar vev (</?) given in (3.2). The fundamental form uj is assumed to have the 
standard local form (2.13). Recall that ip = ipjdz" 1 and x — Xndz 1 A dz 2 are forms on S. 
Substituting in the master equations (3.1) yields 

d 2 ^2 + d 1 ip- 1 -m 2 (z 1 q 1 + z 2 q 2 )x = 

Bix - m 2 (z 1 q x + z 2 q 2 ) ipi = (3.3) 
B 2 x - m 2 (zigi + z 2 q 2 ) ^2 = 

where now x — Xi2- The constants (qi,q 2 ) are the C/(l)j's charges of the fermions that 
belong to a representation (R,qi,q 2 ) of Gs x Z7 (l)i x U(l) 2 . In the following we will 
analyze the different possibilities for the fermions with charges and corresponding curves 
shown in table 1. Notice that the condition d 2 ipi = Biip 2 , implied by the additional 
constraint dip = 0, is automatic by virtue of the last two equations above. 



go, (Ql,Q2) = ( e b°) 

After substituting the charges in (3.3) we obtain the solutions 

^5 = ; x = f(z 2 )e- x ^ 2 ; ^ = -_^ x (3.4) 

where /(z 2 ) is a holomorphic function of £ 2 . The equations (3.3) require the constant Ai 
to satisfy 

Aj = elm 4 (3.5) 

We see that there are solutions localized at z\ — provided that we take the positive 
root Ai = e\m 2 . We then have two zero modes tpi and x which correspond to massless 
fermions of D = 6 massless hypermultiplets living on E . In the presence of magnetic 
fluxes through E tt , chiral four dimensional fermions will appear coming from tpi or /and x 
as dictated by index theorems. 

The characteristic width of the Gaussian wave functions is v — \j[e\m 2 \ We will 
assume that v = l/M 2 , where M* is the F-theory mass scale. For sufficiently large 
compactification radius R this width becomes negligibly small. 

^b, (gi,g2) = (0,e 2 ) 

In this case the solutions of (3.3) turn out to be 

fa = ; x = 9{z l )e~ x ^\ 2 ■ ^- 2 = -^l-x (3.6) 

e2m z 

with g(zi) a holomorphic function of the longitudinal coordinate z\. The constant A2 now 
satisfies 

\\ = ejm 4 (3.7) 
As in the previous situation, having solutions localized at z 2 = requires A2 = e2m 2 . 

gc, (91,92) = (~ei,-e 2 ) 

To treat this curve it is convenient to introduce new variables and fields, and to simplify 
by setting e\ = e 2 = e. Consider then the definitions 

w = z l + z 2 ; ^ = ^(^1 + ^2) 

u = zi - z 2 ; i>u = -(ipi - ^2) (3.8) 



The zero mode equations then become 

2(d w il>w + d u ipu) + em 2 w x = 

du,x + em 2 w ip^ — (3.9) 
duX + em 2 w ipu = 

Now there are localized solutions at w — 0, namely 

^ = ; X = M«)e _AsMa ; ** = \x (3-10) 

em 2 

with h(u) a holomorphic function of the coordinate u along E c , and A 3 = em 2 Notice 
that ipu = implies 

ifo = iJ- 2 = -4= e" A3|w|2 (3.11) 
These agree with results in [33]. 

3.2 Zero modes with fluxes 

We now want to solve the zero mode equations with a background flux, still keeping the 
adjoint vev (ip) given in (3.2). We already know that without flux each curve Ej supports 
localized modes with U(l)i charges given in table 1. The fermions on each curve will now 
feel a total flux T that includes various contributions. There is a bulk U(l) flux F in Gs 
with generator Q (for example, hypercharge or Qb-l)- There are also fluxes along 
the U(l)i inside G p with generators Qi. The total flux can then be written as 

T = FQ + F^Q l + F^Q 2 (3.12) 

The corresponding gauge potentials will be denoted A, A and A^ , with the total potential 
A decomposed like the total flux. We will use conventions in which the covariant derivative 
8a of a field of charge q is defined as 

d A = d-iq A (3.13) 

All field strengths and gauge potentials are taken to be real. 

The fermions x an d 4> have U (1), charges (qi , q 2 ) and transform in some representation 
R of Gs- The bulk flux break G s to T s x U(l) and decomposes into a direct sum of 



irreducible representations that can be labelled by (r,q,qi,q 2 ), where q is the bulk Z7(l) 
charge. The zero mode equations for the charged fermions then become 

(<9 2 - iA 2 )ip 2 + (di - iAi)ipi - m 2 (ziqi + z 2 q 2 ) X = 

{Bi -iAi)x-m 2 {z iqi + z 2 q 2 )*P l = (3.14) 
(B 2 -iA 2 )x-m 2 (z 1 q 1 + z 2 q 2 )^p 2 = 

Clearly, the total gauge potential that appears depends on the charges. It is explicitly 
given by 

A = qA + qi A^ +q 2 A& (3.15) 

The task is to solve the above equations for particular fluxes. 

The 8-dimensional equations of motion further require the vanishing of the (2, 0) and 
(0, 2) components of the field strengths. We will only consider diagonal components T\\ 
and T 22 , even though T\ 2 and T\ 2 are also allowed. Using local coordinates the bulk flux 
takes the form 

F = Fn dz x A dz x + F 22 dz 2 A dz 2 (3.16) 
For the U(l)i fluxes we instead take 

= F 2 { 2 1} dz 2 A dz 2 ■ F (2) = F$ dz x A dz x (3.17) 

The rationale is that, say F^\ is the flux along the curve S a that is defined by Z\ = 
and has coordinate z 2 . 

We will start by analyzing constant field strengths in section 3.2.1. In this case it is 
possible to solve the zero mode equations exactly. We will then study variable fluxes that 
turn out to be necessary to generate corrections to Yukawa couplings [19]. 

3.2.1 Zero modes with constant flux 

In the case of constant field strengths the bulk flux can be written as 

F = 2iM dz 1 A dz x + 2iN dz 2 A dz 2 (3. 18) 

where M and iV are real constants. As explained before, the U(l)i fluxes have components 
only along the curves. They are then given by 

= 2iN {1) dz 2 A dz 2 ■ F (2) = 2iM {2) dz x A dz x (3.19) 



with and some real constants. 

For the gauge potentials we take the following gauge 

A = %M {z\dz\ — z\dz\) + iN (z 2 dz 2 — z 2 dz 2 ) 
A (i) = iN^ (z 2 dz 2 - z 2 dz 2 ) (3.20) 
A® = iM {2 \z x dz x -z x dz x ) 

Notice that the total gauge potential defined in (3.15) can be cast as 

A = iM. [z\dz,\ — Z\dz\) + iN [z 2 dz 2 — z 2 dz 2 ) (3-21) 

where the total flux coefficients are given by 

M = (qM + q 2 M {2) ) ; M = (qN + giiV (1) ) (3.22) 

where q and g« are the bulk and U(l)i charges respectively. In appendix A we give the 
exact solution of the zero mode equations (3.14) with this total gauge potential for the 
three matter curves E , S 6 and S c . Using these results we can then describe the localized 
wave functions at each curve. 

As explained in appendix A, it is convenient to perform a gauge transformation 
A — A + d£l such that A\ = A 2 = 0, and then work with the potential A. We will 
refer to this choice as the holomorphic gauge. The wave functions in this gauge, denoted 
X and tpi, take a simpler form and are better suited to compute gauge invariant quantities 
such as Yukawa couplings. 

In the case of S a we find wave functions 

X = f(z2)e-*W ; & = 5 4 = (3.23) 

where 



M 2 1 M 2 

Ax = —M + e x m 2 Jl + ^— = -M + em 2 + + • • • (3.24) 

e|m 4 2 dm 2 

which reduces to Ai = eim 2 when M. = 0. For future purposes we record the expansion 

of the zero modes to first order in A4, namely 

x = x W{l + M\z 1 \ 2 + ■■■} ; ^ = ^{1- Mv + M\ Zl \ 2 + ■■■} (3.25) 

where v = l/(eim 2 ). Clearly, X ^ = —ipf^ = f(z 2 ) e^ Zl ^ 2 ^ v is the solution for M. = 0. 



Notice that as expected the flux has the effect of deforming the wave function. In the 
holomorphic gauge defined above the wave functions depend on fluxes only through A4. 
Since the matter fields in the curve S a have q 2 = the wave function depends only on 
the flux in the bulk (e.g. from hypercharge in SU(5) or U(1) B -l in SO(10)). Concretely, 
we must replace M. above by 

M a = q a M (3.26) 

where q a is the bulk charge and M comes from the bulk flux. This is relevant later when 
extracting the £7(1) charge dependence of the Yukawa couplings. 

Analogous results are obtained for the E fe matter curve with the obvious replacements 
A4 — > M and ei — > e 2 - In the holomorphic gauge the wave function depends only on 
the bulk flux. This means that Af must be replaced by Af = qbN . 

For the S c curve the wave functions are found to be 

X = h(u + 1 w)e- x ^ 2 e^ ; ^ = \x 5 4 = - (3.27) 
where h(u + jw) is an holomorphic function of its argument and 

£ . M(M - A) 

< = - 3 • (= (A 3 + A) (3 ' 28) 

where A = (M + Af)/2 and A 3 is given in appendix A. In the absence of fluxes one has 
£ = 7 = 0, and A3 = em 2 / y/2, recovering the fluxless result. Note that now it is the linear 
combination £^ + \3ipu which vanishes. On the curve S c the matter fields have U(l)i 
charges (q±, q 2 ) = (— e±, — e 2 ). Hence, M. and M in this case are explicitly given by 

M c = q c M - e 2 M (2) ; N c = q c N - dN^ (3.29) 

We see that the wave function depends on both bulk and brane fluxes. 

3.3 Zero modes with variable fluxes 

In general it is quite complicated to obtain the exact wave functions for non-constant 
field strengths. In [19] an adiabatic hypothesis is assumed whereby the wave functions 
basically follow from those of constant fluxes by replacing the constant coefficients Tq 
by their variable counterparts. An expansion in powers of the z^'s is then performed. 
In this article our approach will be to consider variable fluxes expanded in powers of the 



local coordinates from the beginning, and then solve the differential equations for the zero 
modes. 

We will first expand the fields strengths up to second order in the local coordinates. 
We again turn on only components T\\ and T 22 . Specifically, we take 

T x \ = 2iM + 4i(aizi + o^zi) + U(faz\ + faz 2 ) 

T 22 = 2iM + 4i(a 2 z 2 + a 2 z 2 ) + 6i(fazj + (3 2 z 2 2 ) (3.30) 

where the flux coefficients ctj and fa are complex constants while Ai and M are real. In 
practice the expansion parameter is Zi/R, where R is the compactification radius (see 
section 3.4). We have neglected quadratic terms proportional to (2^) because they do 
not give any new effects concerning Yukawa couplings. The total flux coefficients can be 
split into bulk and curve contributions in analogy to (3.22). 
In our gauge choice the vector potential has components 

Ai = —iMzi — i{(X\z\ + 2a 1 z 1 Zi) — i(faz\ + 3faziZ 2 ) 

A 2 = -iUz 2 - i(a 2 zj + 2a 2 z 2 z 2 ) - %{fi 2 z\ + 3p 2 z 2 z^) (3.31) 

whereas A\ — A\, and A 2 = A 2 . In appendix A we discuss the solutions of the zero mode 
equations (3.14) with this total gauge potential. 

We have not solved the zero mode equations exactly. Instead we found solutions 
in a perturbative expansion in the flux parameters (JA, Af, ccj, We first go to the 
holomorphic gauge with A\ = A 2 = 0, and then iterate to obtain \ — Yl i=$X^\ where 
X* 7 ) is of order / in the flux coefficients. The zeroth order wave function is the fluxless 
solution derived in section 3.1. Once x is determined it is straightforward to deduce the 

For example, in E a , fa = Bix/(eim 2 Zi), and fa = 0. 

The iteration can be carried out to any desired order, but the number of terms will 
clearly be increasingly larger. In appendix A we only display results at most up to second 
order in the flux parameters. Already at first order there is an interesting feature that 
deserves further elaboration. To simplify the argument we set fa — 0. Then, the wave 
functions in the curve S a are found to be 

X = X (0) |l + ^«i-2i+-Mki| 2 + ^i| 2 (ai^i + 2a 1 zi) + ---| 

4 2 1 

fa = ijf ) {l-Mv--va 1 z 1 +M\z l \ 2 + -\z l \ 2 {a l z 1 + 2a 1 z l ) + ---\ (3.32) 



where v — 1/ (eim 2 ) . For constant Tij we have derived the exact solutions whose expansion 
to first order in At agrees with the above results setting ctj = 0. 

One point we wish to make is that in the presence of variable field strength T the 
solution is not merely obtained by adiabatically including the coordinate dependence in 
T. In our case this would correspond to substituting Ai by the effective value 

Meff = M +2(a 1 z 1 + a l z l ) (3.33) 

Indeed, once we replace Ai by At e ff in the solutions (3.25) for constant field strength, 
we reproduce some terms in the expansions (3.32). However, in \ there is an additional 
piece linear in z\ which cannot arise in the adiabatic approximation. In the expansion of 
■01 the term linear in Z\ is expected because in the exact solution there is actually a linear 
term in the constant Ai. 



3.4 Evaluating the fluxes 

Before going to the explicit computation of the Yukawa couplings let us evaluate the size 
of the expected U(l) fluxes in F-theory grand unification schemes. Flux quantization 
demands 

! ' /• m : — / F« = n« ; — f F& = (3.34) 



2tt ./ EC 5 2tt Jx a 2vr Js 



where m, and generically denote flux quanta for the bulk and U(l)i fluxes. On 
the other hand, the GUT gauge coupling constant is given by 

— = Ml [ cuAlu = Vo\(S)M* = R A M^ . (3.35) 

OiG Js 

where R is the overall radius of the manifold S. We then estimate for the fluxes 

F = ■ = 27v^M*ti {1) ; F (2) = 2vr v^M 2 m (2) . (3.36) 

Here we have assumed that the volume of each matter curve is controlled by the overall size 
R, since they are embedded in S. Recall that standard MSSM gauge coupling unification 
gives «g ~ 1/24, for the conventional gauge group normalization TrT 2 = 1/2, with 
generators T in the fundamental of SU(K). Thus, the compactification scale is only 
slightly below the F-theory scale M*. 



Equipped with the above estimates we can characterize more precisely the parametriza- 
tion of the field strengths. For instance, we conclude that the total constant coefficients 
are generically given by 

M = 2n v ^( y qm + q 2 m {2) )M 2 ; M = 2n^/a^(qh + gin (1) )M 2 (3.37) 

where q and are respectively the bulk and U(l)i charges. Similarly, for the total linear 
coefficients we can write 



M 2 \ I ~ 

oti = 2n ^/(^(qotiB + (hot?)—^- ', "2 = 2iry/a^(qa 2B + Qi^)—^ (3.38) 



where olib and 6t£^ are adimensional constants that come respectively from bulk and U(l)i 
fluxes. 

Recall that on the curves S a and the effective wave functions, in the holomorphic 
gauge, depend only on parameters given by bulk quantities. Specifically they are functions 
of 

M 2 

M a = 2n^/a^q a mM 2 ; a la = a a = 2n v fa G ~q a a 1B -^ (3.39) 

M 2 

M b = 27Ty/a^q b nM 2 ; a 2b = a b = 2ny/a^q b a 2B -^- (3.40) 

si 

Other coefficients such as say, M b = 27Cy/^(q b fh + e 2 m ( - 2 ^)M 2 , do not appear in the wave 
functions in the holomorphic gauge A. On the other hand, the parameters for the curve 
S c depend on bulk and U(l)i fluxes according to (see appendix A) 

M c = 2ir^G~(q c m - e 2 m (2) )M 2 ; M c — 2it^/a^(q c n - eih^)M^ (3.41) 
«ic = 2iiyJa G {q c aiB - e 2 a\' )— ; a 2c = 27r^a G {q c a 2B - e x a 2 ) — 

J~C Jri 

In the following we will use A = (M c + M c )/2 instead of M c , and 5 = (a\ c + a 2c )/2 in 
place of a 2c , and we will denote a± c = a c . The decomposition of the quadratic and higher 
order coefficients of T is completely analogous. Observe that gauge invariance imposes 
constraints such as M a + M b + M c = 0. 

Note that the bulk and U(l)i charges, q, q±, and q 2 , depend on the normalization of 
the gauge coupling constants. Consider for example the case of a bulk hypercharge U(l)y 
with q Y integer normalization such that qy(QL, U, D, L, E) = e Y (1, —4, 2, —3, 6). Then, 
TtQy = (12 + 18)e y = 30e y , evaluating at a SU(5) 5-plet. In order to get the standard 




Figure 2: Triple intersection of three matter curves. 

577(5) normalization with TrT 2 = 1/2, a normalization factor ey = l/v60 ~ 0.13 is 
needed. The same exercise for U(1)b~l m SO (10) yields a factor = 1/V24, with 
assignments qB-h(Qh, Qr, L, R) = e B _ L (l, -1, -3, 3). 

The normalization factors for qi^ are found in an analogous way, taking into account 
the enhanced gauge symmetry at each matter curve. Consider for example matter curves 
at which an SU(5) symmetry is enhanced to 50(10) or SU(6). This means that there 
are branchings 50(10) = SU(5) x [/(l)i, 2 or SU(6) = SU(5) x £/(l)i, 2 . One finds 
normalization constants 1 / v^40 and 1 / v^O respectively, with matter fields having charges 
±1. In the case of 5O(10) with matter curves enhancing to Eq or 50(12) one finds 
1 / -\/20 and 1 / respectively. These factors must be taken into account in the explicit 
computation of coupling constants. 

There is an additional constraint on the fluxes in the bulk coming from the BPS 
condition in eq.(2.8) which now reduces to ui AF = 0. In particular, locally this condition 
implies M + N = for constant F. Nevertheless, in what follows we will not impose this 
constraint so that we can keep track of the effect of all flux parameters. 



4 Yukawa couplings 



4.1 Computing Yukawa couplings 

We are interested in evaluating the Yukawa coupling of three chiral fields coming from 
three intersecting matter curves locally described by z\ = 0, z 2 = and z\ + z 2 = in the 
surface S. The piece of the superpotential relevant for Yukawa couplings will be 

W Y = Ml [ Tr (A A A A $) (4.1) 
Js 

where M* is the typical mass scale characteristic of the supergravity limit of F-theory. 
The Yukawa couplings are obtained as overlap integrals over S of the three wave functions 
involved. In principle such a computation requires a knowledge of the wave functions over 
the whole complex surface S. On the other hand, we know that the wave functions are 
peaked around the local curves Z\ = 0, z 2 = and Z\ + z 2 = so that the coupling is 
dominated by the region around the origin, z\ = z 2 = 0, where the three curves meet. 
If this is correct, a local knowledge of the wave functions of the type discussed in the 
previous sections will be sufficient to evaluate the Yukawa couplings. We will thus be 
interested on overlap integrals of the form 1 

Y = Mt \d 2 z x d 2 z 2 ^ 2 ^ 12 (4.2) 
Js 

involving zero modes of the curves z\ = 0, z 2 = 0, and z\ + z 2 = respectively, Zi 
being the local coordinates. Given this structure it is natural to assign the physical, e.g. 
quark/lepton, fields to ipj zero modes and if\ 2 to the Higgs boson and we will assume this 
in what follows. 

We will then take ip\, ip 2 , and ip i2 to be the zero modes localized at the curves S , £5 
and S c respectively. By supersymmetry the wave function of ip± 2 is equal to that of Xc- 
We have seen in the previous sections that, in the holomorphic gauge, the relevant wave 
functions in the presence of variable fluxes take the general form 

fa = -f{z 2 )e-^ m2 \^ 2 G a {z l)Zl] q) 

i>2 = -g(z 1 )e- e ^ 2 ^ 2 G b (z 2 ,z 2 ;q) (4.3) 

<Pi2 = e ^ Gc{zi,z 1 ,z 2 ,z 2 ;q,q 1 ,q 2 ) 
1 See however the note added at the end of the paper. 



where Gi, I = a,b, c, are functions which can be computed to any desired order both 
in fluxes and derivatives of fluxes. Recall that G a and G^ depend only on bulk charges 
whereas G c depends on all charges. In absence of fluxes one simply has G a — G^ — G c — 1. 
Here f(z 2 ) and g(zi) are holomorphic functions. As in [19] we will choose a basis in which 
they are given by fk = (z 2 /R) 3 ~ k and g e = (^i/i?) 3 "^, k,£ = 1,2,3, corresponding to the 
three generations of quarks and leptons. Since there is only one Higgs field we can take 
the corresponding holomorphic function to be a constant. We then have to perform the 
integral 

Y U = Ml [ d 2 z,d 2 z 2 e -«» a (W a +l*l a +£l*+«l a ) f k{z2 ) gfa) G ( Zl , z u z 2 , z 2 ) , (4.4) 

where G = G a GbG c . To simplify the analysis we have set t\ = e 2 = e. 

An important point to remark is that turning on fluxes does not induce mixing in 
the wave functions among different flavors. Indeed, as seen e.g. in eq.(4.3), the flux 
corrections in G a {Gb) do not introduce additional holomorphic dependence on z 2 {z\) which 
would signal generation mixing in the wave functions. This is an important simplification 
because otherwise we would need an additional diagonalization of wave functions in order 
to extract the physical couplings from eq.(4.2). 

The measure in the Yukawa integral can be thought to be 

dn = d 2 z x d 2 z 2 e — 2 (M 2 +NI 2 +^I^I 2 ) (45) 

Clearly, the third exponential, due to the zero mode on E c , breaks the symmetry under 
separate U(l) rotations Z{ — > e *Zi. Instead, there is only invariance under the diagonal 
U(l). This is enough to show that without non-constant fluxes, in which case G can- 
not depend separately on the antiholomorphic variables, the only non-vanishing Yukawa 
coupling is Y33 because f% — g% — 1. Thus, the heaviest third generation of quarks and 
leptons will acquire masses through Y33. 

As pointed out originally in [19], to generate non- vanishing Yukawa couplings for all 
families it is necessary to turn on non-constant background fluxes. To see this it is useful 
to rewrite the measure as 

dfi = J d 2 u d 2 w e - e - 2 (5M 2 +^M 2 ) (4. 6 ) 

where s = \/2 — 1. As before, w = z\ + z 2 and u = z\ — z 2 . In presence of variable fluxes 
the function G can furnish adequate powers of w and u so that the integrand becomes 



invariant under separate phase rotations of w and u. The couplings Yu will thus be 
non-zero and the light generations will gain masses and mixings. We have introduced 
the parameter s in order to study also the case s — 1, which corresponds to ignoring 
the zero mode exponential from E c . In this situation the measure becomes invariant 
under separate U(l) rotations Zi — > e l9l z,i and there will be additional cancellations when 
computing the integrals. 

In performing the integration we will assume that the width of the matter curves is 
determined by the F-theory scale M*, this means v = l/(em 2 ) = 1/M 2 . Consistency of the 
local analysis requires the matter curves to be well localized within S. This amounts to the 
condition v/R 2 <C 1, which is approximately valid for v = l/M 2 because 1/R 2 = y/acM 2 , 
and ^/og — 0.2. In practice we will evaluate the integrals over S, with the above measure 
dfi, by extending \w\ and \u\ to infinite radius. The main contribution to the integrand 
still comes from the region near the origin because the measure is sufficiently peaked. The 
upshot is that in the end all integrals can be done exactly. 

Without varying fluxes there is only one non-vanishing Yukawa coupling Y 33 for the 
third generation which may be explicitly estimated as 

4 0) = Ml [ dW* 2 e- M * (|2l|2+l * 2|2+ > 1+ * 2|2) = • ( 4 - 7 ) 

Js (1 + V2) 

To get the physical Yukawa coupling we really need to work with wave functions nor- 
malized to unity, but to actually normalize our wave functions we would need a global 
knowledge of them over all S. We can however make an estimate by neglecting the effect 
of fluxes and computing the norm of the ipj wave functions from 



Mt 



f e -2M? W ' = * M 2 R 2 (4. 8) 

Js 2 



>s 

Thus, the normalized ipj wave functions are obtained multiplying our wave functions by 
the normalization factor C' 1 ^ 2 . Similarly computed, the normalization for X12, arising in 
the curve E c , is found to be C = C/y/2. We then obtain the normalized third generation 
Yukawa coupling 

Y:a -(i + v ^)Um;W - (1 + V5) "° " 023 (49) 

where we have taken the SU(5) value a G = 1/24. The a^ 4 dependence in eq.(4.9) was 
previously noted in [16]. 



The Y33 Yukawa just computed is given at the unification scale. Taking into account 
QCD renormalization effects down to the weak scale there is an extra factor ~ 3 in the 
case of quarks so that one obtains 

m t ^ 0.23 x 3 x (H u ) = 0.69 x 170 sin/5 ~ 117 GeV (4.10) 

where in the last step we have assumed a large value for tan/3 = (H u /Hd). This is in 
reasonable agreement with experiment, given the uncertainties. A large value for tan/5 is 
required to understand within this scheme the relative smallness of the masses of b-quark 
and r lepton compared to the top. For them one finds 

m b ~ ; m T ~ — 4.11 

tan p 3 tan p 

which gives reasonable agreement for tan f3 ~ 35 — 45. Note that the tau lepton is lighter 
by a factor ~ 3 due to the absence of QCD renormalization. 

There are also subleading contributions to Y33 from flux corrections which appear even 
for constant flux (see appendix B.l). We will eventually neglect all subleading corrections 
so for consistency we will only keep the leading term in Y 33 . When we calculate the rest of 
the Yukawa couplings we will then normalize them relative to the 3rd generation Yukawa 
in eq.(4.7). 

4.2 The case of a constant Xc wave function 

We study first this simple case because it has some interesting features by itself. Further- 
more, a constant wave function is unlocalized and hence could serve to give an idea of 
the results to be expected for Yukawa couplings in which the third particle, presumably 
the Higgs field, lives in the bulk rather than in a localized matter curve. Such type of 
couplings do appear in type IIB and F-theory models in which the base S is not del Pezzo. 
When Xc is a constant, taken equal to one, the Yukawa couplings are determined by 

Yu | X c=i = Mt [ d 2 Zl d 2 z 2 e~ em 2 (l*il 2 +l**l 2 ) f k (z 2 ) g e ( Zl ) G a (z u Zl ; q)G b (z 2 , z 2 ; q) (4.12) 

where q denotes the bulk charges. Substituting the expressions for G a and Gb, which may 
be extracted from the wave functions in appendix A. 2, leads to 

Y M \ Xc= i = 7T 2 4„ 3/ _3 (4.13) 



Hence, the flux-induced distortion of the wave functions does not give rise to any new- 
couplings, only the coupling Y" 33 which is there already for constant fluxes is non- vanishing. 
This is true for any order in the flux expansion. In the next section and in appendix B 
we will provide some examples of the cancellation in the expansion of the Y M . The result 
can also be proven analytically. In fact, notice that in (4.12) the integrals in z\ and z 2 
decouple so that it suffices to show that In — Jd G a vanishes when £ = 1,2. 

The key point is that e~ em2 \ Zl \ 2 G a can be written as — d\F a , as explained in appendix A. 2. 
The function F a can be extracted explicitly, in particular it goes to zero when \zi\ — > oo 
and to 1 when \zi\ — > 0. It is then easy to show that I e = for £ = 1, 2. 

The main conclusion is that in order to get non-trivial fermion mass hierarchies one 
needs all three wave functions to be localized on matter curves. We then proceed to this 
most interesting case of three overlapping localized wave functions. 

4.3 Yukawa matrices 

The physical Yukawas are obtained evaluating the overlapping integral in eq.(4.4) which 
is dominated by the region close to the intersection point. The heavy task is to compute 
the function G by substituting the wave functions found in the previous sections expanded 
in powers both of the flux and derivatives of the flux. In the end each Yukawa coupling 
reduces to a sum of Gaussian integrals that can evaluated analytically. As expected, Y^ 
and Yji are related by an appropriate exchange of flux parameters. 

To begin we have considered the simplest case in which the field strengths are expanded 
only to linear order. This means that we only take into account the first derivative of the 
fluxes (i.e. the oti and 5 parameters) and neglect the effect of higher derivatives. The 
integrals can be determined exactly. For example, the coupling Y" 2 3 is found to be 

v 2 

Y 2 3 = T5 [(s- lfa a + (s 2 - l)a b - 3sa c + s(2s + 3)6] (4.14) 
ori 

where s = y/2 — 1 is the parameter appearing in the measure (4.6). This coupling is 
normalized with respect to Y$ = tt 2 s. This exact expression also shows that when s = 1 
the terms that depend purely on a a and a b completely drop out. In all couplings it 
happens that for s = 1 all pieces involving only parameters of the curves S a and S fe do 
cancel out. This implies that when Xc — 1, the only coupling that survives is Y 33 . 



In appendix B.l we display the leading terms in the expansion in a-fluxes for each 
entry of the Yukawa matrix, normalized with respect to Y^ . Some of the elements Y^ 
have complicated expressions in terms of the flux parameters but the pattern behind can 
be easily understood. Schematically, the couplings turn out to be 

ST < 4 - i5) 

As explained in section 3.4, we have for instance a a = 2iTy/aGq a aiBM 2 / R. Therefore, we 
find v 2 a a /R = 2iraGq a &iB, because v = l/M 2 and M 2 R 2 = l/y/aQ. 

More generally, the corrections to the Yukawa couplings due to first derivatives of the 
fluxes have the general form 

Yij ~ & (2na G (aq + a'q 1 ))^ (2na G (aq + a'q')f- ]) (4.16) 

Here we have simply replaced the atiB and the 6t£^ of section 3.4 by generic constants a 
and a' in order to get an idea of the structure. The ^ are numerical coefficients appearing 
upon integration which are typically in the range 0.1 — 10, as may be seen in appendix 
B.l. The constants q and q' are the bulk and matter curve U(l) charges respectively. 
Recall that for the fields in matter curves S a and S 6 , which include quarks and leptons, 
one has q' — and the corresponding a a and parameters only depend on the bulk 
U(l) charges. This is not the case for the matter curve E c , the parameters a c and 5 do 
depend on the matter curve charges. Note that the normalization of the U(l) charges 
is relevant here. As we explained, for the SU (5) case for integer hypercharge there is a 
normalization factor 1/V60 and the ?7(l)'s on the matter curves containing 10's and 5's 
have normalization 1 / v^40 and 1 / v^60 respectively. 

We have just discussed the general form of each of the terms in the Yukawa couplings 
shown in appendix B.l. To get more accurate results we would need to specify the different 
flux parameters for the three matter curves involved. In particular, we would need a 
precise knowledge of how the U(l) field strengths vary in the vicinity of the intersecting 
points. In principle, given a set of assumptions about the derivatives of fluxes on the 
different matter curves in a concrete model, the formulas in appendix B will allow us to 
compute the different Yukawa couplings. 

It is already quite encouraging that a hierarchical structure of fermion masses seems 
to be built in. Using eq.(4.16) we can further estimate the Yukawa couplings by taking 



into account the normalization of the U(l) charges explained in section 3.4. To this end 
we will write the bulk charges as q = qeB, where q is an integer and is the bulk 
charge normalization. We will similarly write q' = q'e, where e is the normalization of the 
appropriate U(l) on the matter curve, and reabsorb the ratio e/e# into the coefficient 
a'. It is also convenient to introduce cxu(i) — ^B a G, which corresponds to the U(l) fine 
structure constant normalized for integer charges of massless fields. We then conclude 
that the Yukawa matrix has the form 

( ^4 ^3 e 2 v 2 \ 

11 1 
y 6 7] 67] 1 J 



e = 27r v /a c/(1 ) ; r] = ^~a G ~ (4-17) 



In the terminology of [19], we can say that the physical parameter e, which is tied to the 
oti coefficients, controls the flux expansion. The parameter tj is related to powers of the 
width v and the overall radius R that appear in the couplings and controls instead the 
derivative expansion. Taking ac = 1/24 and e# = ey = l/y/60 gives e = 0.16, r] = 0.20 
so that erj ~ 1/31. We then find 

/ "/ \ (3— i) / / ~/ \ (3— j) 

aq + a q V ' ( aq + a q N 



^Hirj (irj (418) 

Therefore, the fermion hierarchies are roughly of the type 

(m 3 : m 2 : mi) ~ (1 : 10~ 3 (ag + a'q') 2 : 10~ 6 (ag + a'q') 4 ) (4.19) 

in qualitative agreement with the observed spectra of quarks and leptons. In the next 
chapter we discuss in slightly more detail to what extent this structure may be successful 
in describing the pattern of quark and lepton masses. 

Let us now see what happens if further terms in the derivative expansion of the fluxes 
were non-negligible. In particular, we have studied the corrections to the Yukawa cou- 
plings when the second derivative flux parameters f3i and p are non-zero. We found that 
the couplings Y 13 and F 2 2 receive leading contributions linear in fa or p. They also have 
quadratic corrections, proportional to the constant coefficients M. and Af of the various 
curves times the fa or p, that are subleading and can be neglected. The leading linear 



terms are 

3 

Y13 = -^[3(s-l) 3 (3 a + 3(s + l)(s 2 -l)(3 b -16s 2 (3 c + 2s(3s 2 + 8s + 6)p] (4.20) 

3 

Y 22 = ^[3(s-l)(s 2 -l)((3 a + (3 b ) + 2s(s 2 -2)p] (4.21) 

where s = \[2 — 1 as before. Here we notice again that when % c = 1 the couplings will 
vanish identically because in this case s = 1 while f3 c = and p = 0. The couplings 
Y12 and Y u have leading corrections typically proportional to ajf3j and f3if3j respectively, 
but the exact expressions are too long to display. In appendix B.2 we show the numeric 
results for the extra leading contributions 

To figure out the size of the corrections due to the second derivative flux parameters 
we will estimate them for the case of the Y n and Y 22 Yukawa couplings. We have 

v v 3 P 27rv 5 ^Wm(bq + b'q')M 2 (h ^«^ (A oo\ 

y 22 ~ W 4 = 2n V^m a G( b( i + b( i) ( 4 - 22 ) 
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Yn ~ -jfr* — ™^ q -^-^ = (2ir) 2 a 2 G a u(1) (bq + b'q') 2 (4.23) 

These terms would contribute to the hierarchy of fermion masses as 

(m 3 : m 2 : mi) ~ (1 : (27r)a Gv /a7^iy : (27r) 2 a 2 G a u(1) ) . (4.24) 

We can evaluate the remaining couplings in the same way (the first order in fluxes Y23 
and Y32 as well). Including the zeroth order Y33 we obtain the structure 
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Since e ~ r], there are hierarchies (1 : e 3 (bq + 5'<f) : e 6 (bq + b'q') 2 ). We see that for 
coefficients of order one, these corrections will generically dominate over the corresponding 
terms in the flux expansion with only first derivatives of fluxes. 

We can go one step beyond and consider also the effect of terms of order three and 
four in the derivative expansion of the fluxes. In this case, to leading order there appear 
contributions to Y 12 , Y 21 and Y" n . The results are given in B.3 where we also explain the 

notation. The additional corrections may be approximated by 

v 5 D 2irv 5 ^/a^(dq + d'q')M 2 „ 

Yn ~ V ( y q -^ = 2n^- ) a 2 G {dq + d'q') (4.26) 
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Y 12 ~ y 21 ~ ^ = v " - = 2n^J^af{cq + c'q') (4.27) 



In this case there are new contributions to the Yukawa couplings of the form 
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(4.28) 



Thus, the first generation Yukawa has corrections of order Y n ~ erf . 

The contributions leading to the terms captured by Y^> are an explicit evaluation of 
the flux expansion of the authors of [19]. On the other hand, their derivative expansion 
would correspond to taking the terms linear in e in F«, Y^ and F (3 ' 4) above. Thus, this 
derivative expansion has the structure 
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(4.29) 
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Note however that, for instance in Y^ 2 \ we also find corrections which do not correspond 
to either of both expansions. 

As a general conclusion, one observes that, for a given Yukawa matrix element, the 
correction due to a higher order term in the derivative expansion will always dominate 
over the flux expansion. 



5 Fermion Yukawa couplings in F-theory GUT's 

In the previous sections we have studied the zero modes of the 8-dimensional quasi- 
topological field theory as well as the computation of the Yukawa couplings using these 
zero modes, without specifying any particular geometry nor identifying the nature of the 
three particles involved in the couplings. We did neither specify the bulk U(l) to be 
considered, which could be hypercharge in SU(5), B-L or other in SO (10). In a GUT 
model, the quark, lepton and Higgs superfields will be localized in matter curves like those 
we have described (see figure 3). We will have Yukawa couplings from a superpotential 
of the form 

W Yuk = Y"Q l WH u + YffQ i D i H d + Y^UWH d (5.1) 

in an obvious notation. In principle the intersection points of the different matter curves 
will be different and, correspondingly the flux parameters 07, (3j, etc., will also be different 



Figure 3: Intersecting matter curves and Yukawa couplings. 



at each intersection. The geometry of each given model may constrain the possibilities 
though. For example, in an SU(5) GUT the left-handed leptons U and the right-handed 
.D-quarks live in the same matter curve. In other settings that need not be the case. 
For example, in a flipped SU(5) setting, D, U and L masses come from independent 
couplings 10 x 10 x 5h, 10 x 5 x 5# and 1 x 5 x 5#. 

Another point to emphasize is that in previous chapters we have made use of the 
possibility of choosing a local basis of holomorphic wave functions of the canonical form 
l,Zi, zf at the intersection point. Note however that in the case of quarks we cannot make 
use of the freedom to choose that basis both at the intersection point leading to [/-quark 
masses and that giving rise to -D-quark masses. On the other hand, if the holomorphic 
basis at both points are very different, one expects very large CKM mixing angles. This 
may be an indication that both points must be quite close in S in order for the basis to be 
aligned to give reduced mixing, as required phenomenologically [19]. That points towards 
further unification into at least E 7 at the F-theory level. 

An important issue in F-theory grand unification is the generation of appropriate 
Yukawa couplings for the [/-quarks. After all, one of the main motivations for going to 
F-theory GUT's instead of perturbative IIB orientifolds is that in the former case these 
couplings are allowed, while they are perturbatively forbidden in type IIB orientifolds. In 
SU(5) the [/-quark Yukawas come from couplings 10* x 10 J x 5#, and if such coupling 




Figure 4: Matter curves intersecting to provide [/-quark Yukawas. 

comes from three distinct matter curves, there can be no diagonal [/-quark couplings. 
This implies that the trace of Y u vanishes which makes impossible a hierarchy of U- 
quark masses. In [15, 16, 19] it was suggested that the matter curve associated to the 
10's in SU (5) could self-pinch as in figure 4-a, allowing for diagonal entries. It was noted 
though [33] that in such configuration the two branches of the wave functions of the 10 are 
independent so that there would be two independent rank one contributions to Y u . This 
would then lead to a rank two Yukawa matrix, with no automatic hierarchical structure. 
In [33] (see also [19]) it was suggested that in fact the two independent branches of the 
wave functions could be identified by some symmetry in the geometry (figure 4-b). In such 
a case the two rank one contributions would be identical and rank one (before the addition 
of flux effects). It was also argued that these symmetries are ubiquitous in F-theory and 
correspond to non-trivial monodromies. Another alternative in order to obtain diagonal 
entries in [/-quark Yukawa couplings was also suggested in [20]. It is more easily described 
in 5*0(10) but it also applies to SU(5). In the context of S'O(IO) the Yukawa couplings 
come from terms 16* x 16- 7 x 10#. If one associates both 16's in the coupling to two matter 
curves S a and and one allows for appropriate B-L flux with opposite restriction on 
the curves, the massless spectrum splits as in figure 4-c. One curve has matter content 
2Q + U, and the other 2U + Q, and this splitting allows for diagonal couplings. In fact, 
as already noted in [20], both matter curves could be local branches of some self-pinched 
matter curve, as in figures 4-a and 4-b. An interesting feature of this possibility is that the 
mixing of the first generation with the other two is expected to be suppressed. In what 



follows we will not specify the particular scheme for the understanding of the [/-quark 
Yukawa couplings. We will just consider that the approximate rank one structure already 
assumed in the previous sections does apply. 

In this section we make a preliminary analysis of the application of our previous 
results to the description of quark and lepton spectra. We would like to see to what 
extent flux distortion may explain the data. Let us first consider for definiteness the case 
of a SU (5) GUT broken down to the SM by fluxes along the hypercharge direction. Let 
us first see whether the flux-induced distortion of wave functions due to first derivatives 
of fluxes is enough to describe the observed structure of quarks and leptons. In order to 
get manageable results we will first assume for simplicity that the fluxes going through 
the third matter curve are approximately constant, i.e. a c = j3 c = 5 = p = (we will also 
denote the subscripts (a, b) as (1,2) hereafter). Under these circumstances the formulas 
in appendix B.l substantially simplify. In particular, the diagonal entries reduce to 

Y 22 ~ eV (0.067(Y^ + Yla\) + 0.11 Y R Y L a x a 2 ) (5.2) 

Yn ~ ey(0A0(Yy i + Yl4) + 0mY R Y L m2(Y^al + Yl4)+0A2YlY^al4) 

where Y R , Yl are the (integer) hypercharges of right-handed fermions (on matter curve 
S ) and left-handed ones (on matter curve E b ). Recall that 

27T 

e = 27rV«y = ^=^c^ ~ 0.16 ; r] = ~ 0.20 (5.3) 

V60 

where we have taken into account the hypercharge normalization factor e Y = l/y/60. 
Note that, as we mentioned before, the wave functions in matter curves S a and E& in the 
holomorphic gauge have no dependence on the U(l)i fluxes, they only depend on the bulk 
fluxes which go in this case along hypercharge. Here a\ )2 are the adimensional constants 
parametrizing the variation of the bulk flux close to the intersection point. Note that in 
principle these parameters may be different for the three different Yukawa couplings, i.e. 
there are a^P' L . To get an idea of the size of the Yukawas let us for the moment assume 
that a(' D,L ~ a 2 ,D,L — a u,D,L . Note that in this case that we neglect the flux variation for 
the Higgs matter curve, the Yukawa matrix is strongly dependent on the hypercharge of 
the quarks and leptons involved in the couplings. Since the maximum value of the quark 
and lepton hypercharges is |y^, aa; j =6,4,2 respectively for leptons and U and D-quarks, 
one expects larger effects for leptons, [/-quarks and D-quarks in that order. 



Inserting the values of the hypercharges for the different particles involved in each 
Yukawa coupling leads to the results in table 2 for the diagonal Yukawas. The couplings 
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0.7(a u ) 2 x 10- 3 
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F^exp) 




(3-4) x 10~ 3 
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5.9 x 10~ 2 
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0.6(a D ) 2 x IO- 3 


3.2(a D ) 4 x 10- 6 


Y D (exp) 




(1 - 3) x to" 2 


(0.6-1.8) x 10~ 3 



Table 2: Hierarchies of fermion masses from the flux expansion. First order in derivatives. 

are normalized to the one of the corresponding third generation particle. We also show 
for comparison experimental results for that hierarchy evaluated at the electroweak scale 
from [46]. The results for the [/-quarks hierarchies are encouraging, for values a u , a L ~ 1 
one can describe reasonably well the observed pattern. For the case of charged leptons the 
mass of the electron is again well described for a L ~ 1. However, the mass of the muon 
would turn out too light unless a L ~ 7.3, which would be quite large and incompatible 
with the a L ~ 1 required for the electron. Thus, the correct numerical description would 
require some further contribution for the muon. Alternatively, it could be that for charged 
leptons neglecting the flux variation coming from the Higgs matter curve is not the correct 
assumption. For the case of the the .D-quark hierarchies one would need slightly large 
values a D ~ 5.8 and a D ~ 4 for Y 22 and Y n respectively. 

Let us now explore what would be the effect of higher order terms in the derivatives 
of the fluxes. If we consider second order in derivatives there are extra corrections which 
may be extracted from appendix B.2. We will again set to zero all flux parameters from 
the curve S c . The corresponding diagonal terms are found to be 

Y 22 ~ erj 2 0.18(y*6i + y L 6 2 ) (5.4) 
Y u ~ eV (-0.23(y^ 2 + y L 2 ^) + 0.29Y L yA5 2 ) 

Again taking b± ,D,L ~ b 2 ' D ' L ~ b u,D,L yields contributions as in table 3. Note that here 
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Table 3: Hierarchies of fermion masses from the flux expansion. Second order in derivatives. 

the Y22 entries have the same structure because in the three cases Yr + Yl — ±3 (we are 
ignoring the overall sign of the contribution which is not relevant for this estimate). As 
expected, the corrections to the Yukawa couplings are always higher than those coming 
from only first derivatives. This is true even for the leptons, which have the highest 
maximal hypercharge and hence get the largest contribution to first order in derivatives. 
In fact, to avoid too large Y^{ L values one rather needs b u ,b L < 1. On the other hand, 
the contribution to Y 2 \ is still too small. The same happens with the D-quarks, one 
would need b D ~ 5.6 to reproduce the observed D-quark mass hierarchies, so that strong 
variation is again required for D-quarks. Terms of order 3 and 4 in flux derivatives could 
also add to the relatively large values of the D-quarks. Equation (4.26) shows that the 
expected contribution is of order Yu ~ erj A {dq + d'q') ~ 2.6c? x 10~ 4 , which reproduces 
the D-quark mass result for a flux parameter d ~ 2.4. Hence, if we do not want to rely 
on relatively large flux parameters the case of D-quarks requires substantial input from 
higher orders in the derivative expansion, up to order four. 

In [19] it was pointed out that the flux expansion to first order in derivatives gives a 
good explanation of the hierarchies observed for leptons and [/-quarks but terms coming 
from the higher derivative flux expansion were needed in order to describe the hierarchies 
for D-quarks. It was also suggested that a possible reason for this different behavior 
could arise from the fact that leptons and [/-quarks have higher maximal hypercharge 
than the D-quarks. We indeed find that the hierarchies for [/-quarks may be quite well 
described by first order flux variations of order one. The resulting electron mass is also of 
the correct order. However, the dependence on hypercharge does not seem to explain the 
different behavior of L and U compared to D fermions. In particular, higher derivative 
terms always generically dominate over the first order terms, even taking into account 



the hypercharge dependence. The milder behavior of the .D-quark hierarchies can be 
understood either by assuming a relatively strong first/second order flux variation (i.e. 
a D ~ 5.8 or b D ~ 5.6) or larger 4th order contributions with d ~ 2.4. The muon has the 
tendency to come out too light which may indicate that neglecting flux variation in the 
Higgs matter curve could perhaps be inappropriate for the leptons and possibly for the 
-D-quark matter curves. 

If we assume that [/-quarks get their Yukawas already at first order in derivatives 
(eq.(4.17)) and on the contrary the _D-quarks need a dominant contribution at order two 
or higher (eq.(4.25) or eq.(4.29), it does not matter for this approximation), we can also 
give an estimate of the CKM mixing matrix [19]. Indeed in this case the respective mass 
squared matrices will be proportional to 
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Then, as in [47], one can estimate the matrices V U,D which diagonalize each of them 
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The CKM matrix, V 
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0.006 0.03 1 j 

which is in reasonable agreement with experiment. This structure is similar to that found 
in [19], although in comparison, in the above formula the separate dependence on the 
hypercharge flux is explicit and the 3rd generation mixing is slightly smaller. 

As a general conclusion, in this simplified scheme in which we have set the flux variation 
in the third curve to zero, one can reproduce the general pattern of quark and lepton 



hierarchies as well as quark mixing, for reasonable choices of flux variation parameters. 
This is particularly the case for the [/-quarks and the electron. Nevertheless, a more 
complete numerical study, not neglecting flux parameters of the Higgs matter curve, may 
be required to get full agreement. The order of magnitude estimates for the CKM matrix 
are on the other hand quite promising. We leave a more detailed phenomenological 
analysis of this framework for future work. 

6 Final comments 

In this paper we have studied the local structure of zero mode wave functions of chiral 
matter fields in F-theory compactifications. We have solved the relevant differential equa- 
tions for the zero modes which were derived from local Higgssing in the world-volume 
effective action of the F-theory 7-branes [15]. These wave functions have a Gaussian pro- 
file centered on the matter curves and become distorted in the presence of U (1) fluxes both 
on the bulk and on the matter curves themselves. In our approach we first write the fluxes 
in a power series of the local coordinates and then make a perturbative expansion of the 
wave functions in powers of the flux coefficients. In this way we obtain expressions which 
may then be applied to compute physical quantities of interest. In this paper we have 
concentrated on the calculation of Yukawa couplings but the wave functions could also 
help to examine other problems. For instance, they could be used to explore the effects 
of closed string fluxes and warping on the effective action, which could prove important 
in relation to compactifications with broken supersymmetry. 

With the wave functions at our disposal we have computed Yukawa couplings by 
performing explicitly the overlap integrals of the three wave functions linked to fermions 
and the Higgs field. By choosing an appropriate gauge, the wave functions of quark 
or lepton generations are shown to depend only on the bulk fluxes but not on the extra 
f/(l)'s associated to the unfolding of the singularities. For example, in the case of a 577(5) 
F-theory GUT broken to the SM by hypercharge flux, the effective distortion of the wave 
function depends on the hypercharge of the specific particle considered. The Yukawa 
integrals can be done analytically and in appendix B we provide the leading terms in the 
flux expansion. One interesting fact we find is that for a constant non-localized Higgs 



wave function, presumably corresponding to a Higgs field living on the bulk of the base 
S, the flux distortion cancels in such a way that the possible Yukawa matrices remain of 
rank one. On the other hand, when the three wave functions are localized, corresponding 
to three intersecting matter curves, a non-constant U(l) flux gives rise naturally to a 
hierarchy of Yukawa couplings as first pointed out in [19]. 

We have applied our findings to the understanding of the observed hierarchies of quark 
and lepton masses and mixings. In a simplified situation in which the flux variation in the 
Higgs matter curve is negligible we obtain explicit compact formulas for Yukawa couplings 
as a function of flux parameters and the charges of the bulk U(l). In a SU(5) setting 
broken to the SM by hypercharge flux, the resulting Yukawa couplings depend on different 
powers of the hypercharge of each quark and lepton. It turns out that reasonable values of 
flux parameters, involving only a first derivative expansion of the fluxes, can account for 
the hierarchical structure of the masses of U -quarks and the electron. The explanation of 
-D-quark hierarchies seems to require larger contributions from the higher order terms in 
the flux derivative expansion. A reasonable semiquantitative understanding of the CKM 
matrix is then obtained somewhat analogous to the results in [19]. 

The natural appearance of hierarchies for masses and mixings looks quite promising. 
However, a full explanation of the data would require a more detailed phenomenological 
analysis. In particular in the numerical estimations we assumed weakly varying fluxes in 
the Higgs matter curve, which needs not necessarily be the case. Furthermore, we also 
took flux variations of the same order for the matter curves corresponding to left- and 
right-handed fermions, which again is suggestive but not generally true. We think that 
our explicit formulas are a good starting point for a more thorough investigation which 
we plan to carry out elsewhere [48]. 

Another interesting topic to address is the origin and structure of neutrino masses, 
which seem to follow a pattern quite distinct from that of quarks. Here the crucial point is 
the nature and origin of the mass of right-handed neutrinos. We think that our results will 
also be useful in this case. More generally, U(l) fluxes may have meaningful implications 
for other physical issues such as supersymmetry breaking. As an example, in [18] it was 
proposed that in F-theory or type IIB orientifolds, local volume modulus dominance of 
supersymmetry breaking gives rise to a very predictive pattern of soft terms consistent 



with radiative electroweak symmetry breaking. It was also pointed out that the presence 
of U(l) fluxes affects in a small but significant way the values of the soft terms and that 
these flux contributions could be needed in fact in order to obtain the proper amount of 
neutralino dark matter. Corrections coming from hypercharge fluxes could also play an 
important role in the detailed understanding of gauge coupling unification [25,35]. It thus 
appears that the distortion caused by fluxes could be indeed important in several physical 
issues in F-theory unification. 
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Note added 

The Yukawa couplings among fields on curves S , T, b and S c arise from the superpo- 
tential term 

W Y = Ml j Tr (A a AA b A $ c ) + cyclic permutations (6.1) 
Js 

where A and $ are chiral superfields given in (2.3). It is enough to focus on 99 terms 
involving two fermions and one scalar. The three families of quark and leptons are taken 
to reside in curves S a and S 6 while the Higgs lives on S c . Then, neglecting an overall 
constant, the coupling is given by 

Yij = [ Wal^yPc ~ ipUriiVc + Mlatfi (6.2) 
Js 

Note that the Yukawa computations in the main text of the paper involve only the contri- 
bution from the first two terms. On the other hand, in a fully symmetric local interaction 
the additional four terms from cyclic permutations should also be included. This has 
been recently addressed in refs. [49] and [50]. In [49] it has been shown that does not 
receive corrections when fluxes are turned on. We wish to stress that there is a delicate 
cancellation among the six contributions in eq.(6.2), each term being in general flux de- 
pendent. This happens independently of whether or not the field strengths satisfy the 
BPS condition lu A F = 0. 

It is instructive to consider the example of constant fluxes. In this case it can be exactly 
shown that each non-trivial term in (6.2) separately gives a flux dependent contribution 
to the third generation coupling Y" 33 , but the full coupling is flux independent. As in 
[49], using our notation, we turn a gauge field along the Q\ and Q2 directions given by 
A = A a Qi + A b Q 2 . Furthermore, we choose 

A a = A b = —iMzidzi — iNz 2 dz 2 + c.c. (6.3) 

Notice that the gauge field acting on S c is A c = —(A a + A b ). The resulting field strength 
satisfies the BPS condition provided M + N = 0. 



The zero modes on each curve follow from the results in appendix A. We find 
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notation. Observe that in the example of constant fluxes these expressions lead to simple 
results such as t/Vi — —v\ a ip a , ip a 2 = 0, and so on, so that only three of the terms in 
(6.2) are not zero. It is straightforward to show that the coupling vanishes except when 
i — j — 3, and that 



Therefore, Y 33 = ir 2 v 2 , independent of fluxes. However, notice that each separate term in 
(6.5) depends on fluxes even if the BPS condition M + N = is satisfied. 

In the example of [49], in which the BPS condition u A F = is satisfied, it also 
happens that the flux effects on the couplings only cancel when all terms in (6.5) are 
included. On the other hand, in the setup of this article, in which tv A F = is not 
enforced, nonetheless it can be checked that when all terms in (6.5) are added only the 
coupling Y" 33 survives and is flux independent. In [50] the sum of all contributions to the 
couplings has also been taken into account. 

In [49] it was proved that the cancellation of flux effects in the full coupling follows 
from an exact residue formula. For a pedestrian derivation of this formula we start from 
(6.2) and manipulate the integrand to write it as a sum of total derivatives. To this 
purpose, following [49], we write the zero modes ipij and ipi, which satisfy the last two 




(6.5) 



where the ipi are given in (6.4). Evaluation of the Gaussian integral yields 
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equations in (A. 2), as ipij = dj£i, together with 
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The functions hi are holomorphic and correspond to ipj . An elementary calculation, 
dropping family indices to simplify, then shows that 



Y = d 2 Z 1 d 2 Z 2 {dl [(lp a 2(fib ~ VWa)£c] - d- 2 [{ipalVb ~ AWa)£c}} 
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The integrals in the first line can be evaluated by parts, and then the boundary terms are 
seen to vanish because the zero modes ip a and ifb are localized. In the second line, inte- 
grating by parts twice, using £ a = v(tp a — h a )/zi, similarly for £&, and invoking localization, 
gives the final residue formula Y ~ Res ( ha z hb z hc J [49]. 



The computation of Yukawa couplings just described is purely local. If the symmetry 
among the cyclic permutations in eq.(6.2) still remains after a global completion of the 
theory, only one generation acquires a Yukawa coupling. In this case the observed hier- 
archy of fermion masses cannot be generated just by turning on magnetic fluxes, some 
additional ingredient, e.g. non-perturbative effects, should also be at work to produce 
these mass hierarchies. 




A Fluxed zero modes and wave functions 



In this appendix we study the solutions of the zero mode equations (3.14) both for constant 
and variable field strengths. We will explicitly consider the curves E a and E c . As in the 
fluxless case, the results for E a and E& are completely similar, but the curve E c must be 
treated separately. 

We find it convenient to rewrite the total gauge potential as 

A = A + dVl (A.l) 

in such a way that A\ = A 2 = 0. We can then work in this 'holomorphic' gauge where the 
potential is just A and the corresponding fermions are denoted x an d i>- The advantage 
is that the equations reduce to 

(<9 2 - iA 2 )ip 2 + (di - - m 2 (zi<?i + z 2 q 2 ) X = 

dix - n?{z x q\ + z 2 q 2 ) i>\ = (A.2) 
d 2 x - vrP^qx + z 2 q 2 ) i) 2 = 

and the gauge fields do not appear in the last two equations. The further constraint 
8x>P = becomes dip = and is automatically verified on account of the last two equations 
above. 

The solutions for the original flux are recovered by performing a gauge transformation, 
namely 

X = e in X ; ^i=e"Vi ; ^ 2 = e l ^ 2 (A.3) 

To compute Yukawa couplings it suffices to work with the hatted fields because the cou- 
plings are gauge invariant. 

A.l Constant flux 

From the total gauge potential given in (3.21) it follows that the transformed potential 
and gauge function are 

A = -2iMz x dzx - 2iHz 2 dz 2 ; Q = i(M\zr\ 2 + Af\z 2 \ 2 ) (A.4) 

We then need to find the solutions of (A.2) when A\ = —2iMz\ and A 2 = —2iMz 2 . The 
charges (qi,q 2 ) that must also be specified depend on the curve. 



go, (Ql,Q2) = ( e b°) 

Notice that in this case M = qM and M = qN + eiA^ 1 ), where (M, AT) come from the 
bulk flux and N^ from the flux along the curve. As in the fluxless case, we find that 
tjj 2 = 0, which then implies d^x — 0- We make the Ansatz 

X = f{z 2 )e-^\ 2 (A.5) 

The equation d\x = zi?Pi then fixes 

k = — (A.6) 
There is still an equation that requires Ai to satisfy 

\\ + 2M\ X - efm 4 = (A.7) 
To have localized solutions we choose the root 



M 2 

A, = -M + r.nr + ( A - 8 ) 

which reduces to Ai = e\m 2 when Ai = 0. Inserting in (A.5) and (A.6) gives the solutions 
found in [19] in a different gauge. 

gc, (gi,g2) = (~ei,-e 2 ) 

In the fluxless case we saw that to solve the equations it is convenient to set e\ = e 2 = e, 
and to use the variables w — (z\ + z 2 ) and u — (z± — z 2 ), together with the redefined 
fermions tjj w = (ipi + ip 2 )/2, and ipa = (ipi - ip 2 )/2. 

The gauge potential is still formally given by (A. 4) but now Ai = (qM — e 2 M^) and 
AS = (qN - eiAf (1) ). In the new variables the non-vanishing components of A are 

Aw = -iAw - i(M - A)u 

A u = -iAu - i(M - A)w (A.9) 

where A = (Ai +M)/2. In the gauge A, the zero mode equations imply that the t[j 
fermions neatly depend on x as 

iv = — \- du,x ; fa = — \— daX (A.10) 
em z w em z w 



In turn \ can be determined from the remaining equation 

(d w - iA w )iJu, + (d u - iA u )iJu + ^em 2 wx = (A- 11 ) 
To solve we make the Ansatz 

X = h(u + jw) e- A3|tu|2 e*™ (A.12) 
where h(u + jw) is a holomorphic function of its argument. It then follows 

i>w = — ; i>u = 2 x (A.13) 

em z em z 
Substituting in (A. 11) determines the unknown constants. We find 

Finally, A 3 is a positive root of the cubic equation 

A 3 (A 3 + M)(X 3 + 2A-M) - \e 2 m\\ z + A) = (A.15) 

When A4. = Af = we recover the fluxless solution with £ = 7 = 0, and A 3 = em 2 / y/2. In 
the special cases A = (Af = —M.) and A = M. (Af = M) the cubic becomes quadratic 
and the positive root is easily identified. 

A. 2 Variable flux 

We consider the quadratic flux given in (3.30). The corresponding transformed potential 
and the gauge function turn out to be 

Ai = -2iMz x - %{pL]z\ + 2a 1 z 1 z 1 ) - 2i(fi x z\ + 3(3^1) 

A 2 = -2iAfz 2 - 2i(a 2 z 2 2 + 2a 2 z 2 z 2 ) - 2i(fo4 + ?>{5 2 z 2 z\) (A.16) 

n = i\z x \ 2 [M + %{axz x + a x z x ) + (fi x z\ + A~z\)\ 

+ %\z 2 ^ [Af + (a 2 z 2 + a 2 z 2 ) + (f3 2 z 2 2 + (3 2 z 2 2 )] 

As described below for particular curves, we have only been able to obtain zero mode 
solutions in a perturbative expansion in the flux parameters. 

So, (51,52) = (ei,0) 



As in the constant flux case we find ip 2 — which implies d 2 x — 0. On the other hand, 
ipl = j^dix, where v = \ je\m 2 . There is still an equation 



(di — — eim 2 zi x — 



(A.17) 



with A\ given in (A. 16). We have found a solution x — X]/=o^ ( ' / ' ) > wnere X^ is °f order 



/ in the flux coefficients. There is a corresponding expansion for ip\ with ip^ 



Zl 



diX [I) - 



The zeroth order solutions are the fluxless ones presented in section 3.1. They are 



X (0) = f(z2) e 



-eim 2 |2i| 2 



$°) = -f( Z2 ) e -^?W? 



(A.18) 



The expansion of x to second order turns out to be 



X = X (0) |l + v 3 H 23 + v 2 H 22 + v(H 21 + H u ) + #io + \h 2 + ■ ■ ■ j 



(A.19) 



where i> = 1/exm 2 is the volume defined before. The auxiliary functions are given by 



H 2 \ 



= Z\Z\ 



M + -(aizi + 2a lZl ) + ^fe 2 + Sfrz 2 ) 



4 3 ^ 2 



(A.20) 
(A.21) 



l -M 2 + \M(2a x z x - a lZl ) - {\a\ + %M^)z\ + {- A Mh + ^a 2 )z 2 



16 <_> 4 _ — _o 1 _ „ o _ 2 — o 15 _q 4 ^ 7)9 — 4 

—aiP 1 z 1 + -QJi/Ji^ - -ai/Ji^i^i + —a 1 fj 1 z 1 z 1 - yPi% + gPi z i 



(A.22) 



#22 = ^1 



4. . ,4 9 3 . . _ . 16 . 9 15., , 

y o 4 o 

2 , - _ . 

to 



#23 = 



to 



(A.23) 
(A.24) 



It can be checked that when a.\ — (3± — 0, the results match those of section A.l to second 
order in M. 

The expansion of the wave function ipi needed to compute Yukawa couplings follows 
from tpi = —&iX- We obtain 



ik = i>T\ 1 + v 6 H* z + v 2 K 22 - v(M + H* u + K 2l ) + H 10 + -Hf + ■ ■ ■ } (A.2o) 



1 



with the additional definitions 

K22 



16 .4 2 9 A - . 16 - 2 15-2 3 

+ + -MPJz! + -g-aiM + y/W 



— 2i^i [2a 1 p 1 z 1 - oixfixzx) 
15 



(A.26) 



^21 — ^1^1 



+ ^( 10q;iZi + I3«izi) + (Ja 2 + - A Mfix)z\ + (|a 2 + ^-M(3i)zl 



.4 
9' 



9 
4' 



,4 
3 ( 



11 

T 



i^l |2 - 1 14 — ^ 2- _ 2 ,„ 1 2 2-2 

+ -|Q!i] Zi^i + — a^z-^zi + — a x p\Z 1 z 1 + 3|pi| 2^ 
3 5 5 

+ ^ol^zI + + \p 2 A + j^ 2 ** 



(A.27) 



E c> (gi,g2) = (-ei,-e 2 ) 

We need to solve the zero mode equations (A. 10) and (A. 11). The gauge potential com- 
ponents A w and A u can be easily found changing to coordinates w = (z,\ + z%) and 
u — {z\ — z 2 ) starting from (A. 16). As before we define A = (M + M)/2. In analogy we 
also introduce 

5 = ^(«i + a 2 ) ; P = l -{f3 l +f3 2 ) (A.28) 

and the corresponding 5 = 5* and p = p*. 

To iterate we begin with the zeroth order solutions presented in section 3.1, taking 
h — 1. They are 



£(0) = e -em 2 M 2 /V2 



# 0) 



,-em 2 \w\ 2 /V2 



(A.29) 



To higher orders we will only report the wave function \ that enters in Yukawa couplings. 
To first order we find 



V2 



D u + D 



10 



(A.30) 



with functional coefficients given by 



D w = w 



(M - A)u + - Aw + - (oil - 5) wu + - (an - 5) (wu + uw) + -5 (w 2 + 3m 2 ) 
2 2 2 6 



+- 5 (ww + uu) + - — p) (w 3 + -uw 2 ) + - (Pi — p) (w 2 u + 2uww + 3u 2 u) 

1 3 

+ — p (w 3 + 6-um 2 ) + — p (w 2 w + Awuu + 2-uru 2 ) 
16 16 



= w 



4 9 



(A.31) 
(A.32) 



We have also computed the second order correction to \. We refrain from presenting it 
because it involves too many terms. 



B Yukawa couplings 

The purpose of this appendix is to provide the explicit expressions for the Yukawa cou- 
plings Yij obtained upon performing the overlap integral of the localized wave functions on 
the curves Ej, I — a,b, c. The procedure is to determine the integrand z%~ 1 z\~ 3 G a G b G c , 
where the Gi are the corrections of the wave functions due to fluxes that were derived in 
appendix A. The integral with measure (4.6) is evaluated assuming that the size of the 
compact manifold S is much larger than the width v of the matter curves. The piece of 
the integrand that can contribute is a sum of terms |w| 2m |w| 2n and the integral is easily 
computed. 

B.l Flux expansion, first order in derivatives 

To proceed systematically, to begin we consider the field strengths expanded up to linear 
order. In this way we obtain Yukawa couplings depending only on the parameters olj 
and S that characterize the first derivative of the total flux acting on the matter curves 
Ej. We focus on the leading terms for each entry. There are corrections proportional 
to powers of the zeroth order coefficients, M. and A/", times powers of the oci, which are 
always subleading, i.e. higher order in e. We have normalized with respect to the zeroth 
order third generation Yukawa coupling = n 2 s, where s = \pl — 1. The results are 
as follows: 

v 2 

Y 23 = — [0.11 a a - 0.27 a b - 0.41 a c + 0.52 5] (B.l) 
R 

Y 22 = -^[0.067(a 2 + a^) + 0.11 a a a b - 0.02 a a a c 

+ 0.17 c? c + 0.02 a b a c + 1.11 5 2 - 5(0.33a c + 0A5a b + 0A2a a )] (B.2) 
v 4 

y 13 = — [ - 0.03a 2 - 0.16a 2 ; - 0.32a 2 - 0.04a a at, + 0.02a a a c 
R 2 

- 0.30a fe a c + 0.15a a <5- 1.835 2 + 1. 21a c 5 + 0.97 a b 5] (B.3) 



v 6 

Y 12 = — [ 0.03a 3 + a 2 (0.02a fe + 0.01a c - 0.125) 

+ a a (0.03a 2 + 0.01a fe a c - 0.03a 2 - 0.21a b 5 - 1.08a c 5 + 1.595 2 ) 

- 4.325 3 + 4.5a c <5 2 + l.67a b 5 2 - 0.91a 2 5 - 1.5a b a c S 

+ 0.12a 2 b 5 + 0.20al + 0.08a b a 2 c -0.06a 2 b a c -0.07a 3 b ] (B.4) 

Y n = ?L[o.i(ai + a$) + a 3 (0.09a 6 - 0.15a c - 0.535) + a 2 (3.105 2 - 0.15a c 5 

- 1.03a b 5 + 0.15a 2 - 0.04a 6 a c + 0.12a 2 ) + a a (-1.05 3 - 18.96a c <5 2 + 4.785 2 a 6 ) 
+ <5a a (8.40a 2 - 0.43a 6 a c - 1.11a 2 ) + 0.15a a a 3 + 0.22a a a 6 a 2 

+ 0.04a a a 2 a c + 0.09ajja a + 2.855 4 + 39.31a c 5 3 - 4.07a 6 <5 3 - 21.11a** 2 

- 16.50a fe a c <5 2 + 3.4a 2 <5 2 + 1.46a 3 5 + 9.33a b a 2 J - 0.45a 2 a c 5 

- 0.83a 3 5- 0.36a 4 - 0.15a 6 a 3 + 0.15a 2 a 2 + 0.15a^a c ] (B.5) 

The couplings satisfy the property 

Y i:j (a a , a b , a c , 5) = Yji(a b , a a , 25 - a c , 5) (B.6) 

Then, the Y^ for % > j can be easily found from the above results. 

We want to stress that just as Y 23 given in (4.14), all couplings can be computed 
exactly. The results are given numerically only for ease of presentation. For example, Y22 
is found to be 

Y22 = ^[s(s-l)(s 2 -l) 2 (a 2 a + a 2 ) + ^(s-l) 2 (3s 2 + l)a a a b 

+ s(s-l)(3s-l)(a a -a f ,)a c -^v / 2s(3s 2 + 4v / 2s-4)a 2 (B.7) 

1 3 
+ -s(16s 3 - lSv^s 2 - 40s + SGv 7 ^ 2 + -s(3v / 2s 2 + 8s - 4V2)Sa c 

+ ^s(12s 3 - 15s 2 + 4s - 9)5a a + ^s(12s 3 + 3s 2 - 20s - 3)5a b ] 

where s = \[2 — 1 is the parameter in the measure (4.6). 

For completeness we also provide the expansion of Y33 to first order order in fluxes, 
namely 

Y 33 = 1 + v [ X - (s - 1) (M a + U b ) + s A ] (B.8) 

with s = \/2 — 1. This is the simplest example showing that the corrections vanish when 
Xc = 1 which implies s = 1 and A = 0. 



B.2 Flux expansion, second order in derivatives 

To second oder in derivatives of the fluxes there are further contributions to the Yukawas 
with leading terms as follows: 

v 3 

Yi3 = ^[~ 0.07& - 0.44ft - 0.34ft + l.Olp] (B.9) 

Y22 = ^ [0.18(ft + ft) ~ 0-57p] (B.10) 

v 5 - - - 

Y 12 = — [ft(-0.16a a - 0.04a b + 0.04a c + 0.135) + ft(0.16a a + 0.39a 6 + 0.02a c - 0.615) 

+ ft(0.91a a + 1.3a b + 0.36a c - 4.135) + 5(-0.81a a - 2.00a 6 - 1.7a c + 6.295] (B.ll) 

+ ft(2.03p - 1.74ft) - 5.49ft 2 + 10.99ftp - 7.20p 2 ] (B.12) 

The coupling F31 follows from Y"i 3 by exchanging ft <-> ft and ft <-> (2p — ft). A similar 
remark applies to Y 2 \. 

B.3 Flux expansion, third and fourth orders in derivatives 

In section 4.3 we also discuss the effects of third and fourth order derivatives in the fluxes. 
The modified wave functions needed to calculate the couplings are obtained as explained 
in A. 2 but with new terms in the gauge potential because now the components of the 
total field strength have the additional pieces 

Tlt a = Hd4 + C t zf) + m(D lZ f + Dfzt) (B.13) 

To third order the effective flux acting on the £/ is characterized by parameters Cj and 
A c as discussed in section 3.4. The notation at fourth order is analogous. The couplings 
that receive new corrections are 

5v 3 

Yu = ^[(s-l)(s 2 -l) 2 (D a + D b )+2s(s 4 -3s 2 + 3)A d ] (B.14) 

Y 12 = T ^[Q(s-l) 2 (s 2 -l)C a + Q(s 2 -l) 2 C b -15s(s 2 -2)C d 

+ s(12s 3 + 15s 2 -20s-30)A c ] (B.15) 

Y21 = ^[6(s 2 -l) 2 C a + 6(s-l) 2 (s 2 -l)C b +15s(s 2 -2)C d 

+ s(12s 3 -15s 2 -20s + 30)A c ] (B.16) 



where s = V2 — 1. Observe again that these couplings vanish when 
evaluation gives 

Y n = ^[-0.50(D a + D b )+2.Q0A d )] 

Y 12 = ^[-017C a + 0.41(7 6 + 1.13(7 c -1.44A c ] 

Y 21 = ^[-017C b + 0.4ia-l-13C c + 0.83A c ] 
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